Unexpected Behavior of the Local Compressibihty Near the B=0 Metal-Insulator 

Transition 
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We have measured the local electronic compressibility of a two-dimensional hole gas as it crosses the 
B=0 Metal-Insulator Transition. In the metallic phase, the compressibility follows the mean-field 
Hartree-Fock (HF) theory and is found to be spatially homogeneous. In the insulating phase it de- 
viates by more than an order of magnitude from the HF predictions and is spatially inhomogeneous. 
The crossover density between the two types of behavior, agrees quantitatively with the transport 
critical density, suggesting that the system undergoes a thermodynamic change at the transition. 



PACS numbers: 71.30.+h, 73.40.-c, 05.70.Ce 

The recent discovery of a Two-Dimensional (2D) 
Metal-Insulator Transition (MIT) Q has raised the fun- 
damental question concerning the existence of metallic 
2D systems mQ. In contrast to the scaling theory of 
localization |g|], which predicts that only an insulating 
phase can exist in 2D, there is compelling evidence for 
metallic-like behavior in a growing number of 2D systems 
[Q,^. To date, the vast majority of this evidence comes 
from transport measurements. Clearly, if a true zero- 
temperature phase transition exists, it may reveal itself 
also in the thermodynamic properties of the 2D gas, such 
as the electronic compressibility, k~^ = ri?-^ (where /i 
is the chemical potential and n is the carrier density). 
For example, a crossover to a gapped insulator would re- 
sult in a vanishing compressibility at the transition [p|, 
whereas a crossover to an Anderson insulator would in- 
volve a continuous evolution of the compressibility across 
the transition. 

The compressibility or t^ of a system reflects how 
its chemical potential varies with density. For non- 
interacting electrons it simply amounts to the single- 
particle Density of States (DOS), which in 2D is den- 
sity independent (^ = 1^)- This picture, however, 
changes drastically when interactions are included. Ex- 
change and correlation effects weaken the repulsion be- 
tween the electrons, thereby reducing the energy cost, 
thus leading to negative and singular corrections to j^. 
Within the Hatree-Fock (HF) theory, which includes both 
the DOS and exchange terms, one gets: 
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with the compressibility becoming negative at low 
enough densities. Measurements of the macroscopic com- 
pressibility of 2D electron/hole gases, in the metallic 
regime, have indeed confirmed this behavior Rfl- 

In this work, we have expanded the study of /i(?i) into 
the MIT regime. Our measurements utilize several Sin- 
gle Electron Transistors (SETs), situated directly above 
a Two-Dimensional Hole Gas (2DHG). This technique al- 
lows us to determine the local behavior of /i(ri) as well as 



its spatial variations. Simultaneous macroscopic trans- 
port measurements were conducted to ensure a precise 
determination of the MIT critical density in the same 
sample. Our measurements indicate a clear thermody- 
namic change in /i(n) at the MIT. In addition, we find 
that the behavior oi fi{n) in the metallic phase follows the 
HF model and is spatially homogeneous. The insulating 
phase, on the other hand, is found to be spatially inho- 
mogeneous. In this regime, -^ deviates by more than 
an order of magnitude from the predictions of the HF 
model. 

The samples are inverted GaAs/AlGaAs heterostruc- 
tures p. A layer of AlGaAs, 2000 A thick, separates 
the 2DHG from a p+ layer that serves as a back gate, 
and is separately contacted. A GaAs spacer, 500A thick, 
separates the 2DHG from the doped GaAs layer situ- 
ated just under the surface (see inset in Fig. 1). The 
mobility of the 2DHGs ranges from 7 • 10^ cni^/Vsec to 
1.2 • 10^ cm'^/V sec at n = 5 • 10^" cto^^ and T = 4.2 K. 
Biasing the back gate enables us to vary the density of 
holes by two orders of magnitude, from 2 • 10^^ cm~^ to 
2 • lO^CTO^^. Our samples show a clear MIT. A typical 
plot of the resistivity, p, vs. density at various tempera- 
tures is shown in Fig. 1. A clear crossing point between 
metallic and insulating behavior is observed. The critical 
density and resistivity, determined from the low tempera- 
ture crossing point, is TT-c ~ 2.1-10^° cm~^ andpc ~ 0.9-^. 
Similar samples were previously used to study the trans- 
port properties of the MIT in GaAs |g| . 

A local measurement of /i requires the ability to mea- 
sure electrostatic potentials with high sensitivity and 
good spatial resolution. It has been previously shown 
that both can be achieved using an SET 0-0. We 
have therefore deposited several aluminum SETs on top 
of an etched Hall bar mesa (inset in Fig. 1). This 
configuration allows us to study simultaneously the lo- 
cal thermodynamic behavior and the macroscopic trans- 
port properties and, thereby, trace possible correlation 
between them. Our technique relies on the fact that /i 
of the 2DHG changes as its density is varied. At equi- 
librium, the Fermi energy is constant across the sam- 



pie, and therefore, a change in /i induces a change in 
the electrostatic potential, which is directly measured by 
the SET. The size of the SET and its distance from the 
2DHG determine its spatial resolution, estimated to be 
0.1 X 0.5 fim^ . The measured voltage sensitivity of the 
SET is 10 iiV. Similar technique has been previously em- 
ployed to image the local compressibility of a 2DEG in 



the quantum Hall regime 1 1 1 12] . 

To demonstrate the strength of our technique, we first 
employ it to study the magnetic field (B) dependence of 
fi deep in the metaUic regime (n = 1.96 • 10 



11 



cm 



see 



Fig. 2). Theoretically, the formation of well-separated 
Landau levels at high magnetic fields is expected to in- 
duce sharp saw-tooth like oscillations in fi. The sharp 
drops are at integer filling factors {v), and the slopes 
obey Ih = (j + ^) ^ ± 5*MB, where m*, g* are the ef- 
fective mass and g-factor, j is the Landau level number, 
and the ± relate to the two spin directions. We find 
good quantitative agreement between this model and the 
measurement using a single value for m* = 0.34 • rue and 
g* — 2.7. The measured value of m* is also in good 
agreement with the known effective mass of light holes 
in GaAs, to* = 0.38 ■ TOe IQ. The position of the jumps 
allows us to determine the local density of the 2DHG un- 



der the SET: m = 1.82 ■ 10" 



cm 



which deviates by 



only 7% from the density measured macroscopically by 
transport techniques. This small deviation demonstrates 
that the area beneath the SET is hardly disturbed by its 
presence. 

The main focus of this work is the density dependence 
of /i across the B — MIT. Both the theory and previous 
macroscopic measurements [R0] show that at low densi- 
ties fi should increase monotonically as n is decreased 
(Eq. 1). Unexpectedly, our local measurements demon- 
strate quite a different behavior. Several typical traces of 
fi vs. n are shown in Fig. 3. Fig. 3a shows the behavior 
in the metallic regime and Fig. 3b shows the behavior 
across the MIT and in the insulating regime. Instead 
of the expected monotonic dependence, /i exhibits a rich 
structure of oscillations. All the oscillations, including 
the fine-structure seen on the left side of Fig. 3b, are 
completely reproducible and do not depend on the mea- 
surement rate or sweeping direction. Two distinct types 
of oscillations are observed: In the metallic regime (Fig. 
3a) we observe long saw-tooth oscillations with a typical 



period in density of (1 — 2) • 10 



10, 



". Superimposed on 



them and starting in close proximity to the MIT, a new 
set of rapid oscillations emerges (Fig. 3b). Their typical 
period is an order of magnitude smaller: (1 — 2)-10^ cto^^, 
and their amplitude grows continuously from their point 
of appearance to lower densities. Similar behavior has 
been measured on seven different SETs placed on five 
separate Hall bars from two different wafers. Although 
the specific pattern of oscillations varies between these 
Hall bars and changes also after thermal cycling, they all 
have the same general characteristics. The qualitative 



change in the oscillation pattern of /x(n) near the critical 
density suggests that the system experiences a thermo- 
dynamic change at the MIT. 

Out of the two types of oscillations, the easier ones 
to explain are these on the metallic side. In Fig. 3a 
we compare the measured oscillations with the expected 
monotonic behavior of the HF model. In contrast to the 
model, the measured /j, has a density independent average 
which suggests that some kind of screening mechanism is 
present in the system. While the negative slopes of fj,{n) 
(black symbols) follow the HF theory, there are appar- 
ent additional drops (some of which are extremely sharp) 
between them (gray symbols). This saw-tooth profile 
is reminiscent of the behavior of the chemical potential 
of a quantum dot as a function of its density |14| and 
suggests the existence of discrete charging events in the 
dopant layer, situated between the 2DHG and the SET. 
Thus, the measured /x of the 2DHG varies undisturbed 
along the negative slopes, until a certain bias is built be- 
tween the 2DHG and the SET that makes the charging 
of an intermediate localized state energetically favorable. 
This causes a screening charge to pop-in which result in 
a sharp drop in the measured electrostatic potential, af- 
ter which 11 continues to vary undisturbed until the next 
screening event occurs. Because screening occurs only at 
discrete points we can reconstruct the underlying fj.(n) 
from the unscreened segments in the measurement. We 
do this in Fig. 4, by collecting the slopes of these seg- 
ments from the saw-tooth profiles of five different SETs 
placed on three separate Hall bars. Each single point in 
this graph represent the slope of a well-defined segment, 
like the ones emphasized in Fig 3a. In the metallic side 
all the data collapses onto a single curve (the dashed line 
in Fig. 4 is a guide to the eye). The prediction of the HF 
model is depicted by the solid line in Fig. 4. The finite 
width of the 2DHG adds a positive contribution to the 
compressibility as shown by the dotted line. It is evident 
that both theories describe the data to within a factor 
of two over more than an order of magnitude in density, 
throughout the metallic regime. The fact that different 
Hall bars and different SETs (namely different locations 
within the same Hall bar) produce the same dependence, 
suggests that the metallic phase is homogeneous in space. 

The sudden appearance of rapid oscillations as the sys- 
tem crosses to the insulating side already signifies that 
something dramatic happens at the transition. Assum- 
ing that this new set of oscillations is caused by the same 
mechanism, namely screening by traps, we can proceed 
and extract its slopes and add them to the same plot of 
Sn/Sn (see Fig. 4). Unlike in the metallic phase where 
the system clearly has a negative Sfi/Sn, namely negative 
compressibility, in the insulating phase we a-priory do not 
know the sign of the compressibility. We, therefore, plot 
in the left side of Fig. 4 both negative and positive slopes, 
all of which deviate considerably from the expected n~^' ^ 
power law. The deviation becomes greater than an or- 



der of magnitude at our lowest density ||I5[ . We see that 
6^/5n starts deviating from the theoretical curve in close 
proximity to the transport-measured critical density, sig- 
nifying a change in the screening properties of the 2DHG 
at the transition. Another intriguing result is the non- 
universal behavior of the slopes on the insulating side. It 
should be emphasized that each slope is determined from 
a complete segment in the /i (n) curve (see inset to Fig. 
3b) and is, therefore, determined very accurately. The 
fluctuations in the slopes, seen even in a measurement 
done using a single SET, are completely reproducible and 
suggest that mesoscopic efii'ects are present. Furthermore, 
the average behavior of Sfj,/Sn in this fluctuating regime 
is seen to be position dependent. An example of that is 
shown in the inset in Fig. 4 where we plot Sfi/Sn mea- 
sured by two separate SETs on the same Hall bar. Such 
dependence on position indicates that once the system 
crosses into the insulating phase it ceases to be spatially 
homogeneous. 

The charge traps have a clear signature in our mea- 
surements, emphasizing their important role in the ther- 
modynamic ground state of the system. It was recently 
suggested by Altshuler and Maslov ||l^ that a gas that is 
in electro-chemical equilibrium with charge traps might 
show metallic-like behavior of the resistivity. This be- 
havior is a result of a temperature-dependent scattering 
mechanism |17 evoked by the coupling to the trap sys- 
tem. Although it is clear from our measurements that 
charging of the traps takes place as the density of holes is 
varied, we are currently unable to determine experimen- 
tally whether this coupling is responsible for the metallic 
behavior seen in our samples. 

The large values of 6fi/6n in the insulating phase is 
intriguing. Currently there are no theoretical predictions 
for large negative 6n/6n^ which in turn leads to a sharp 
rise in fj, as the density is reduced. However, we would 
like to dwell on several scenarios in which the measured 
dfi/dn becomes positive at the transition and increasingly 
large as the density is further reduced. The flrst scenario 
assumes the formation of voids in the 2DHG. Such voids 
are unable to screen the back gate voltage and, hence, 
would result in an added large and positive contributions 
to Sfi/Sn. In this scenario, the size and shape of the 
voids, as function of decreasing density, are responsible 
for the large and seemingly random set of positive slopes 
observed. A second scenario involves the formation of 
isolated puddles in the 2DHG iQ . These puddles screen 
the back gate discretely and, hence, would produce large 
and positive slopes of /i (n) . It should be noted that 
both scenarios suggest a percolativc like transition MJ20| 
and are in accord with our observation that the 2DHG is 
inhomogeneous in the insulating phase. 

In conclusion, we have measured the dependence of 
/i(n) across the MIT. Our measurements clearly indicate 
a thermodynamic change in the screening properties of 
the 2DHG at the transition. We find the metallic phase 



to be spatially homogenous and to behave according to 
the predictions of the HE model. The insulating phase, 
on the other hand, deviates significantly from the HE 
predictions and is spatially inhomogeneous. 
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FIG. 1. Density dependence of the resistivity: Traces from 
top to bottom on the right correspond to T=600, 550, 500, 
450, 400, 300 and 250 mK . The low-T crossing point is at 
Uc = 2.1 ■ 10^° cm-'^. Inset: An SEM micrograph of the SET, 
placed on top of a schematic structure along with the mea- 
surement circuit. 



FIG. 2. Measured B-dependence of ^. The solid line is the 
predicted Landau fan with m* — 0.34 ■ rUe and g* = 2.7. 
Oscillations in n are observable up to i^ = 7 . The local 
density extracted from the magnetic field values at integer i/ 
is ni — 1.82 • lO'^^ cm~^ which deviates only by 7% from the 
transport-measured macroscopic density, n = 1.96T0^^ cm^^. 



FIG. 3. (a) fi{n) in the metallic regime. The HF predic- 
tions (y) are depicted by the solid line. The negative-slopes 
are high lighted (dark symbols) to demonstrate their resem- 
blance to the HF model, (b) /i(n) across the MIT and in the 
insulating region. Inset: A closer look at the data in the in- 
sulating regime. Each slope is composed of many data points 
allowing for an accurate determination of the slopes. 



FIG. 4. Sfi/Sn collected from five SETs on three differ- 
ent hall bars from two different wafers. In the insulating 
regime, both negative and positive slopes are shown (closed 
and open symbols respectively). Each point corresponds to a 
well-defined segment in the /.t(n) trace. The point marked by 
an arrow corresponds to the marked segment in Fig. 3a. Inset: 
Results from two SETs on the same device (•, A) demonstrat- 
ing the spatial dependence of Sfi/Sn on the insulating side. 
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